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ABSTRACT 

The  dynamics  of  a  proposed  marine  nuclear  propulsion  plant, 
under  the  condition  of  sudden  changes  in  load,  are  studied 
by  both  an  analog  and  digital  simulation.   Two  inherent  reac- 
tor characteristics,  internal  reactivity  feedback  and  the 
existance   of  delayed  neutron  emitters ,  which  effect  power 
plant  stability  and  controllability  are  considered  in  detail. 

Comparison  of  the  simulation  results  with  practical  data 
indicates  that  the  simulation  represented  the  dynamics  of  a 
marine  nuclear  power  plant  with  a  sufficient  degree  of  accuracy 
to  be  extremely  useful  in  the  study  of  marine  nuclear  power 
plants . 

Evaluation  of  the  simulation  results  indicate   that  the 
proposed  power  plant  is  inherently  stable  when  operating  in 
the  normal  region.   The  desirability  of  control  schemes  are 
discussed  and  a  control  scheme  utilizing  a  constant  average 
coolant  temperature  program  is  implemented.  This  external 
control  scheme  significantly  improved  the  response  of  the 
system. 
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I.   INTRODUCTION 

Nuclear  reactor  power  plants  for  use  in  marine  application 
differ  from  stationary  power  plants  in  that  they  must  respond 
quickly  to  power  demands  which  may  vary  throughout  a  wide 
range  in  a  relatively  short  period  of  time.   There  are  maneu- 
vering situations  where  the  shipboard  power  plant  may  be 
required  to  respond  to  load  changes  of  up  to  50  percent  in  a 
few  seconds.   In  addition  to  the  operational  differences,  the 
mobile  nature  of  the  nuclear  power  ship  coupled  with  the  size 
and  weight  limitations  which  must  be  imposed  upon  the  reactor 
shielding  require  that  safety  be  a  primary  consideration  in 
design  and  control  of  the  marine  power  plant. 

There  are  several  references  [Harrer  1963  and  Glasstone 
and  M.  C.  Edlund  1952]  that  treat  the  nuclear  reactor  and 
nuclear  reactor  control  from  the  standpoint  of  the  reactor 
operating  at  constant  output  power.   Schultz  [Ref.  3]  discusses 
the  problems  of  reactor  control  from  the  standpoint  of  a 
reactor  operating  with  varying  power  demands  as  well  as  when 
operating  under  steady  power  out  conditions.   King  -[-Ref-; — 4-}" 
considers  the  problems  associated  with  nuclear  power  systems 
largely  from  the  position  of  a  mechanical  engineer. 

This  thesis  investigates  a  marine  nuclear  power  plant  and 
is  intended  as  an  initial  look  into  the  dynamics  of  this  type 
of  propulsion  system.   The  system  dynamic  equations  are 
eveloped  and  the  resulting  fundamental  balance  equations  are 
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simplified  by  invoking  reasonable  first  approximations. 
Initially  a  power  plant  without  an  external  controller  is 
studied.   After  the  initial  investigation  a  control  scheme  is 
formulated  and  implemented.   The  plant  response  with  controller 
is  compared  with  the  plant  without  controller. 

The  power  plant,  represented  in  block  diagram  form  in 
Figure  (1-1) ,  utilizes  a  pressurized  water  reactor  fueled  with 
low  enriched  uranium  dioxide,  heat  exchanger  for  the  conver- 
sion of  reactor  power  to  steam  and  a  turbine  which  is  repre- 
sented by  the  changing  load.   The  reactor  and  one  side  of  the 
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Fig.  (1-1) .   BLOCK  DIAGRAM  OF  BASIC  NUCLEAR 
POWER  PLANT 


heat  exchanger  represent  the  primary  loop  while  the  other  side 

of  the  heat  exchanger  and  turbine  represent  the  secondary  loop 

in  the  system. 
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The  system,  although  hypothetical,  is  similar  to  the 
plant  on  board  N.  S.  Savannah.   Many  of  the  initial  (before 
load  changes)  temperatures  and  coefficients  closely  approxi- 


pj® 


mate  those  of  Savannah  f-5-7 — 6-7 — 7i~.   In  order  to  show  the  valid- 
ity of  the  mathematical  model  used  in  this  study,  a  qualitative 
comparison  between  the  response  of  the  simulation  and  both 
those  of  Savannah  and  a  proposed  marine  power  plant  is  used 
[8]. 

The  organization  of  this  thesis  follows  the  format  as 
described  below:   Chapter  (2)  contains  a  development  of  the 
mathematical  model  used  for  the  evaluation  of  system  perfor- 
mance; Chapter  (3)  contains  a  general  discussion  of  two 
inherent  reactor  characteristics,  delayed  neutrons  and  internal 
reactivity  feedback,  that  affect  the  stability  of  the  power 
plant;  Chapter  (4)  presents  an  analysis  of  the  analog  and 
digital  computer  simulations  of  the  power  plant  without 
external  control  and  also  contains  a  comparison  of  the  results 
of  the  two  simulations;  Chapter  (5)  discusses  some  general 
aspects  of  nuclear  power  plant  control  and  shows  the  specific 
control  scheme  used  in  this  evaluation;  Chapter  (6)  is  a 
discussion  of  the  conclusions  arrived  at  by  this  study  and 
also  discusses  areas  of  further  interest.   Appendices  A 
through  D  supplement  these  chapters. 
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II.   MATHEMATICAL  MODEL 

The  mathematical  model  for  the  plant  under  consideration 
is  developed  by  obtaining  the  model  for  each  system  component 

individually.   The  individual  models  together  with  the  values 

\ 
of  the  specific  system  parameters  are  summarized  at  the  end  of 

this  chapter.   During  the  course  of  the  development ,  it  is 

necessary  to  make  several  simplifying  assumptions.   These 

assumptions  introduce  some  inaccuracies  in  the  system  model. 

"However,  the  errors  involved  in  these  assumptions  are 

usually  less  than  the  amount  of  uncertainty  in  the  engineering 

/ 
value  of  the  coefficients  used."  [3]   The  generalized  equa- 
tions given  below  are  listed  in  [Ref.  9], 

A.   REACTOR  KINETICS 

The  point  model  kinetic  equations  of  a  chain  reacting 

nuclear  pile  have  been  developed  in  the  literature  [10,  11] 

and  shall  be  used  for  the  starting  point  in  the  development 

of  the  plant  model  to  be  studied. 

6 

E 

i=l 


-(^1  ■*    *    h*i  w-w 


Bi 

ci  =  —  n  -  XL   ci  (2"2) 


Equations  (2-1)  and  (2-2)  are  state  equations  describing  a 
nuclear  reactor  operating  at  some  steady  state  power  condition 
The  variables  n,  c,  and  6k  represent  neutron  density,  density 
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of  delayed  neutron  emitters  and  excess  reactivity.   The  other 
terms  are  constant  coefficients  described  ©a>-grage"^. 

The  indexed  terms  in  equation  (2-1)  and  (2-2)  indicate  a 
multi-group  representation  of  the  delayed  neutron  emitters. 
This  six  group  representation  is  valid  for  a  detailed  study 
of  reactor  dynamics,  however,  in  the  scope  of  this  study  the 
reactor  is  considered  as  a  power  source  which  is  investigated 
from  the  terminal  point  of  view.   Therefore,  without  loss  of 
generality,  and  for  the  convenience  and  simplicity  (i.e., 
equipment  limitation  in  the  analog  system) ,  the  kinetic 
equations  are  represented  by  a  one  group  delayed  neutron 
approximation.   Equations  (2-1)  and  (2-2)  thus  become: 


6k     3n    ,  /0  _, 

n  =  —  n-ij—  +Ac  (2-3) 


c  =  £*  -  Ac  (2-4) 


The  magnitudes  of  the  variables  n  and  c  are  far  greater 
than  the  magnitudes  of  the  other  system  variables.   In  order 
to  facilitate  the  use  of  the  analog  computer,  equation  (2-3) 
and  (2-4)  will  be  normalized  by  dividing  through  by  their 
steady  state  values  n(0)  and  c(0). 


"/"<°>  -(r1)  Hw-(r)   hw+hw  <2-5' 


c        _     /  3   \      n  Ac  o-^ 

cloT  "    ll  )  3ToT  "  cToT  c      } 
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Next  defining: 


N  =  -5-*-  ;  D  =  -4rr  ;  and  K  =  ~ 
n(o)        cCO)  3 


noting  that  in  steady  state: 

N  =  D  =  1,   N  =  D  =  K=0;  and    n  (0 )  =  X    c(0) 


Equations  (2-5)  and  (2-6)  thus  become: 


N  =  |  (KN  -  N  +  D)  (2-7) 


D  =  X(N  -  D)  (2-8) 

The  term  K  above  is  the  effective  net  reactivity  existing 
within  the  reactor.   In  steady  state  there  is  a  constant 
neutron  population  and  therefore  K  must  be  zero.   To  establish 
K  at  this  zero  value  the  reactivity  attributed  to  rod  position 
must  be  off  set  by  the  inherent  internal  negative  feedback 
(internal  feedback  is  discussed  in  some  detail  in  Chapter  3) . 
Thus  the  algebraic  equation  for  K  is: 


K  =  KR  +  ap  Tp  +  ac  TAV  (2-9) 


Tp  is  the  average  fuel  temperature  and  T    is  the  average 

coolant  temperature.   The  coefficients  a_  and  a   are  coeffi- 
^  F      c 

cients  of  reactivity  feedback.   It  should  be  noted  that  a 
and  a   are  negative  .numbers . 
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B.   HEAT  TRANSFER  EQUATIONS  FROM  FUEL  AT  T   TO  COOLANT 
AT  T^ 

For  a  given  change  in  temperature  the  fuel  will  absorb 
heat  in  the  amount  of 


qf  =  mf  cf  CT2  -  T  ) 


T2  >  Tx 


(2-10) 


Where  qf  is  the  heat  absorbed  by  the  fuel,  mf  is  the  mass  of 
the  fuel,  cf  is  the  specific  heat  of  the  fuel  and  T2  -  T      is 
the  change  in  fuel  temperature.   The  time  rate  of  change  of 
the  heat   absorbed  in  this  manner  is: 


2t  mf  Cf  (T2  "  V 


(2-11) 


Consider  one  reactor  fuel  element  of  Figure  (2-1) .   The 
temperature  profile  of  the  fuel  varies  along  the  axial  length 
Z  and  depends  upon : 

1)  the  heat  input  from  the  reactor. 

2)  the  coolant  flow  rate,  which  is  constant. 


f 


Eh 

2 

W 

J 

2 

w 

W 

D 

& 

* 

Fig.  (2-1) .   TEMPERATURE  PROFILE  OF 
SINGLE  FUEL  ELEMENT. 


The  first  of  several  assumptions  must  now  be  made.   Assume 
that  the  temperature  profile  does  not  change  shape  or  form 
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during  transient  conditions  and  that  in  the  small  signal 
analysis  cf  and  mf  are  constant.   Then  at  any  axial  position 
Z  the  rate  of  heat  absorption  by  the  fuel  is  dependent  upon 
the  time  derivative  of  fuel  temperature.   Further  considering 
changes  from  steady  state  where  T,  is  a  constant,  it  then 
follows : 


qf  =  mf  cf  TF  (2-12) 


The  general  heat  transfer  equation  is  [12] 


q  ■  k,  V2  T  (2-13) 


axial  (z  direction)  flow  can  be  neglected  in  comparison  to 
radial  flow.   For  a  homogeneous  right  cylinder  fuel  element 
equation  (2-13)  becomes: 

employing  the  power  balance  for  each  fuel  element  where  Power 
=  time  rate  of  change  of  heat 


Power  generated  \   /   Power    \     /     Power 
by  the  fission    J  =  I   absorbed   ]  +     transferred 
process       /   \by  the  fuel/    \  to  the  coolant 


Considering  the  case  where  k...  is  independent  of  space 
variables,  the  (differential)  power  balance  equation  across  a 
differential  radius  r  is: 

P  -  mf  cf  Tp  t  ^   2^  +  |  »  (2-15) 
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or  rearranging 


f   m 


f  °t 


m 


2 (   2^T    1  2T\ 

■  o       {  2t2    r  2rj 


(2-16) 


Since  cf  and  k.  are  both,  functions  of  temperature,  the  above 
equation  is  a  non-linear  partial  differential  equation. 
Following  the  assumption  that  the  temperature  profile  does 
not  shape  or  form  and  that  2T/2r  is  a  constant  as  shown  in 
Figure  (2-2) 


TEMP 


COOLANT 

2T    AT 
2r    Ar 


T  -T 

F   AV 

r2-rx 


22T 
2r 


=  0 


RADIUS 


TAV  = 


Tu+T 
H   c 


Fig.  (2-2) .   TEMPERATURE  PROFILE  ACROSS 
FUEL  ELEMENT  AND  CLADDING. 


thus  equation  (2-16)  becomes 


T   =  — 

F   mf  cf 


mf  cf 


T   -  T 
F    AVI 

r2  -  r± 


(2-17) 


The  power  generated  is  a  function  of  the  neutron  population 
[1] .   Thus: 


--—-  =  f  (N)  =  A,  N 

m.p  c,  1 
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kl  -1 

The  term  has  the  units  of  sec    and  can 

mf  cf  rCr2  -  r^ 

be  related  to  a  time  constant  t..   Therefore  equation  (2-17) 


becomes : 


*f  =  Al  "  "  7^  (TF  "  TAV»  .   (2"18) 


C.   HEAT  BALANCE  FOR  COOLANT  IN  THE  CORE 

Under  the  same  assumptions  as  stated  above  the  rate  of  heat 
absorbed  by  the  coolant  as  it  passes  through  the  fuel  elements 
is  : 

*c  ■  mo  cc  at  (th  "  V  (2"19) 

Where  q   is  the  rate  of  heat  transferred  to  the  coolant,  m 
^c  c 

is  the  mass  flow  rate  of  the  coolant,  c   is  the  specific  heat 
of  the  coolant,  T   is  the  hot  leg  temperature  and  T^,   is  the 
cold  leg  temperature.   Employing  a  power  balance  as  above 

/  Power  transferred  \  _  /  Power  absorbed  \   /Power  transferred 
I   to  the  coolant    /    \  by  the  coolant  j       \      to  the  boiler 


k   /  T   -  T   \  m   c 

_±        _t av  \  T   +  _c — c  (T   _  T  j       (2-20) 

r   \  r^  -  r.,   /     c   c   AV      t     H     C 


Where  x   is  the  time  for  a  unit  volume  of  coolant  to  flow 
o 


through  the  reactor. 


k,  (T   -  T  )       . 
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.       kl 

again -—7 r can  be   related   to   a   time   constant  t_. 

r(r0    -   r, )    rn      c  2 

2  lcc 


Thus 


*AV  =   TZ    'TAV  "  V    "   7-    (TH   "   TC»  <2"22> 

2  o 


Rearranging 


T  T  2T 

*AV  "   xf  "   if*  +   T^  <2-23' 

2  2  o 


where 


o         2 

T       +    2t0 
o  2 


D.   HEAT  TRANSPORT  DELAYS 

There  are  time  delays  associated  with  the  transfer  of  heat 
from  one  portion  of  the  plant  to  another.   In  developing  the 
model  for  the  time  delays  the  following  assumptions  must  be 
made : 

1)  there  is  no  heat  loss  in  the  piping  between  the  reactor 
and  the  boiler. 

2)  there  is  no  mixing  of  fluids  in  the  system  piping. 
1.   Time  Delay  From  Reactor  to  Boiler 

The  boiler  inlet  temperature ,    T   ,  is  a  delayed  version 
of  the  reactor  outlet  temperature  Tu. 

n 


TBI  "  TH  (t  "  T3'  (2"24) 
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Where  T    is  the  boiler  inlet  temperature  and  x   is  the 
transport  time  between  reactor  and  boiler.   Rearraning 


TH  "  TBI  Ct  +  T3»  (2"25) 


Expanding  in  a  Taylor  Series 


TH  =  TBI  +  T3  TBI  +  •  •  •  (2"26) 


Solving  for   T_T 


*BI  "  TJ  (TBI  "  2  TAV  +  V  (2"27) 

2 .   Time  Delay  From  Boiler  to  Reactor 
Following  the  above  reasoning 

Tc  ■  77  (Tc  "  2  tb  +  TBI»  (2"28) 

4 

Where  T_,  is  the  average  temperature  in  the  boiler. 

E.   HEAT  BALANCE  FOR  HEAT  EXCHANGER  (BOILER) 

The  development  of  the  model  of  a  counter  flow  heat 
exchanger  in  which  there  is  a  change  of  phase  (water  to  steam) 
in  one  flow  path  is  a  difficult  task.   It  is  not  the  intention 
of  this  chapter  to  present  a  rigourous  derivation  of  the  equa- 
tions used  in  the  system  model,  but  rather  it  is  to  give 
insight  into  the  physical  relationships  from  which  these 
equations  were  obtained.   Therefore  a  more  heuristic  argument 
will  be  used  in  the  development  of  the  steam  generator  model. 
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TEMPERATURE 


Fig.  (2-3) 


DISTANCE 

TEMPERATURE  PROFILE  FOR  COUNTER 
FLOW  HEAT  EXCHANGER. 


There  are  several  additional  assumptions  which  must  now  be 
made : 

1)  There  is  no  heat  lost  or  gained  with  the  surroundings. 

2)  The  specific  heat  of  the  individual  fluids  remains 
constant . 

3)  The  overall  coefficient  of  heat  transfer  (k2)  is 
uniform  through  the  heat  exchanges. 

4)  From  Figure  (2-3)  the  mean  temperature  difference 

throughout  the  heat  exchanges   T   =  T_  -  T„ 
3  3      m     B     S 

1.   Primary-Side 

Consider  a  small  section  of  heat  exchanger  having  a 
transfer  area  dA.   Then  the  incremental  heat  transfer  is: 


dQ  =  k_  dA  AT   =  c   m   dT 
2      n    c   c   c 


(2-29) 


integrating  yields 


Q  =  k2  A[TB  -  Ts]  =  cc  mc  [TBI  -  TB0] 


(2-30) 
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Since  k-,  A,  c  ,  and  m  are  all  constant,  equation  (2-30) 


c 
reduces  to 


2(TB  -  Ts)  =  E(TBI  -  TBQ)  (2-31) 


Where  T   is  the  steam  temperature  and  A  is  the  area  of  the 
heat  exchanger. 

2 .   Secondary-Side 

Energy  balance  in  the  secondary-side  consist  of  the 
amount  of  power  transferred  to  the  secondary,  k,A(TR  -  Tq) , 
equal  to  the  power  absorbed  within  the  boiler  plus  the  power 
which  is  transferred  to  the  turbine.   To  bring  this  relation 
into  a  linear,  ordinary  differential  equation  some  additional 
assumptions  are  required.   These  are: 

1)  The  thermal  capacity  of  the  boiler  metal,  steam 
and  coolant  can  be  lumped. 

2)  The  power  demand  on  the  boiler  can  be  measured  as 
a  fraction  of  normal  load.  ty  will  represent  this 
condition  with  \p   =  1  representing  the  normal  load. 

3)  The  heat  exchanged  in  the  boiler  occurs  at  a  single 
point  within  the  heat  exchanger. 

Therefore  the  flow  of  heat  delivered  to  the  turbine  may  be 
represented  as  the  product  of  the  throttle  opening,  represented 

by  K_ ,  and  the  steam  pressure.   Since  the  steam  pressure  is 

2PC 
upon  the  steam  temperature,  the  term   j^r-       To  can  ^e  used 

S 
to  represent  steam  pressure.   Thus  the  heat  balance  for  the 

secondary  side  becomes: 

2P 
klA(TB  -  Ts)  =  (mc  cc  +  ms  cs)  Ts  +  KT  *   ^-   Tg    (2-32) 
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Where  ni   and  c   are  the  mass  flow  rate  of  the  steam  and  the 
s      s 

specific  heat  of  the  steam  respectfully.   Rearranging 

klA  KT   *  2PS 

S    "   m      c      +   m      c       L    B  SJ         Tin      c +ln      c)  2T_      S 

cess  cess  S 


(2-33) 


Which  reduces  to 


*S  -  -  h    C"TB  +  TS  +  *1   *  2ff  )  (2'34) 

D  S 


F.   FINAL  MODEL  WITH  SPECIFIC  SYSTEM  PARAMETERS 

The  development  of  the  specific  equations  for  the  system 
is  carried  out  in  detail  in  Appendix  A,  they  are  summarized 
here  for  convenience. 

1.  Reactor  Kinetics 

N  =  6.4  (KN-N+D)  (2-35) 

D  =   .1  (N-D)  (3-36) 

K  =  KR  -  1.57  x  10"2  TF  -  3.14  x  10~2  TAV       (3-37) 

where  N(0)=1,  D(0)=1  and  K  (0)=24 . 7275 . 

K 

2 .  Heat  Transfer  From  Fuel  to  Coolant 

Tp  -  -  |  (Tp  -  TAV)  +  30.0  N  (2-38) 

where  T„(0)=565°F. 

r 

3.  Heat  Balance  for  Coolant 


TAV  =  10  TF  "  50  TAy  +  40  Tc  (2-39) 


where  T   (0)=505°F. 
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4.   Time  Delays 
'-   *BI  =  "  ?  (TBI  "  2  TAV  +  V  <2-40' 


where  TDX (0)=520°F. 

oJL 

*c  =  "  I  CTc  -  W3  -  I  V  (2"41) 


where  T  (0)=490°F. 

5 .   Heat  Balance  for  Steam  Generator  Secondary-Side 

Ts  =  -  j    (-  |  TBI  +  (0.667  +  0.0665  i\>)    Tg)      (2-42) 

where  TS(0)=565°F. 

G.   METHOD  OF  SOLUTION 

Because  of  the  multiplication  of  the  time  dependent 
variables  K  and  N  in  equation  (2-7)  and  because  K  is  a  function 
of  N  the  entire  set  of  dependent  equations  which  describe  the 
plant  dynamics  is  non-linear.   Since  (in  general)  non-linear 
differential  equations  do  not  lend  themselves  readily  to 
analytical  solutions,  computer  simulation  will  be  used  to 
solve  these  equations. 
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III.   INHERENT  REACTOR  CHARACTERISTICS 

There  are  several  important  characteristics  of  nuclear 
power  plants  which  greatly  effect  the  stability  and  control- 
lability of  the  system.   The  two  most  significant  factors 
relating  to  the  system  performance  are  the  existance  of 
delayed  neutron  emmiters  and  the  negative  internal  reactivity 
feedback  associated  with  temperature  changes  within  the 
system.   Both  of  these  phenomena  are  discussed  below  in  some 
detail. 

A.   PROMPT  AND  DELAYED  NEUTRONS 

At  the  instant  of  fission  there  are  neutrons  (prompt) 
which  are  produced  instantaneously.   Fortunately,  for  control 
purposes,  a  small  portion  (£)  of  all  neutrons  produced  during 
the  fission  process  are  delayed  in  time  [1,  2].   These  delayed 
neutrons  are  the  result  of  the  decay  of  fission  products.   As 
noted  previously  these  delayed  neutrons  are  grouped  in  six 
distinct  groups  with  different  delay  constants  (A) .   Analysis 
of  the  effects  which  delayed  neutrons  exhibit  on  the  reactor 
kinetics  can  best  be  shown  by  the  evaluation  of  the  kinetics 
equations  (2-3)  and  (2-4)  with  and  without  the  delayed 
neutrons  considered. 

The  analysis  of  the  prompt  and  delayed  neutron  response 
which  follows  is  a  study  of  the  reactor  operating  at  a  power 
level  low  enough  so  that  feedback  effects  (discussed  later 
in  this  chapter)  are  negligible.   An  analagous  situation  in 


28 


a  power  plant  would  be  in  the  start-up  power  range  of  the 
reactor. 

In  this  evaluation  the  input  to  the  reactor  is  in  the 
form  of  a  perturbation  to  the  steady  state  reactivity  k.   The 
input  was  limited  to  a  maximum  value  (8)  of  (0.0064),  for 
when  6k  exceeds  6   the  reactor  is  in  a  response  condition 
described  as  "prompt  critical."   "Prompt  critical"  refers  to 
the  fact  that  the  nuclear  fission  chain  reaction  can  be  main- 
tained by  means  of  prompt  neutrons  alone.   If  this  condition 
occurs  neutron  density,  and  hence  power  increase  rapidly 
from  the  instant  of  input,  making  the  reactor  difficult  to 
control  and  may  cause  what  is  referred  to  as  a  start  up 
accident.   The  "prompt  critcial"  operation  is  avoided  in 
practice. 

1.   Kinetic  Equations  Assuming  No  Delayed  Neutrons 

Under  the  assumption  of  no  delayed  neutrons,  equations 
(2-3)  and  (2-4) ,  which  describe  the  reactor  kinetics  in  steady 
state  with  a  one  group  delayed  neutron  approximation,  reduce 
to  the  single  equation: 

A  =  i^  (3-D 

Solving  for  n 

6k  . 
n  =  n(0)  e  l  (3-2) 

Dividing  by  n(0),  substituting  the  system  parameters  found  in 
Appendix  A,  and  assuming  6k  =  0.00  3  equation  (3-2)  now  becomes 


N  =  e 


3t  (3-3) 
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2 .  Reactor  Kinetic  Equations  with  a  One  Group  Delayed 
Neutron  Emitter 

"*  Under  the  one  group  delayed  emitter  assumption, 

equations  (2-3)  and  (2-4)  are  unchanged.   With  the  substitution 

of  system  parameters  from  Appendix  A  and  from  the  solution  of 

these  equations  carried  out  in  Ref .  [11] ,  the  solution  of 

these  equations  for  N  is : 

N  =  1.88  e+0-0884  fc  -  0.88  e"3"4  t.  (3-4) 

3.  Comparison  of  Results 

The  effect  of  the  delayed  neutrons  is  readily  seen 
by  the  comparison  of  equations  (3-3)  and  (3-4) .   With  no 
delayed  neutrons  the  reactor  neutron  population  has  increased 
by  a  factor  of  8/100  in  3  seconds,  however  when  delayed 
neutrons  are  present  the  reactor  neutron  population  has  only 
increased  by  a  factor  of  2.5  in  3  seconds.   Plots  of  N  versus 
time  for  both  of  these  situations  are  shown  in  Figures  (3-1) 
and  (3-2) .   As  time  passes  (i.e.,  time  greater  than  10  seconds) 
the  first  term  of  equation  (3-4)  becomes  dominant  and  the 
response  will  follow  a  positive  exponential  increase.   However, 
in  the  time  span  of  interest  (a  small  time  after  the  distur- 
bance) the  affect  of  the  delayed  neutrons  is  to  greatly  reduce 
the  rate  of  change  of  neutron  flux.   It  is  the  natural  occur- 
ance  of  these  delayed  neutrons  which  makes  it  possible  to 
operate  a  nuclear  reactor  with  the  degree  of  safety  required. 
Even  as  small  a  fraction  as  0.6  4  percent  of  the  total  neutron 
population  acting  as  delayed  neutrons  has  the  effect  of  making 
the  entire  reactor  control  problem  greatly  simplified  [2,  3]. 
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It  is  considered  pertinent  to  re-emphasize  that  the 
above .discussion  of  prompt  and  delayed  neutrons  was  conducted 
for  a  reactor  in  the  startup  range.   The  response  of  a 
reactor  operating  in  the  normal  power  range  (where  reactivity 
feedback  is  strong)  is  significantly  different  than  that  of 
the  reactor  operating  at  low  power.   This  difference  will  be 
demonstrated  later  in  this  chapter. 

It  should  be  pointed  out  however,  that  the  delayed 
neutron  phenomenon  also  has  an  undesirable  facet.   While  it  is 
readily  obvious  that  it  restricts  the  growth  of  the  neutron 
population  for  positive  6k  inputs,  it  also  restrains  the 
reduction  rate  of  neutron  population  which  would  result  from 
a  negative  6k  input  as  the  power  demand  is  reduced  during 
normal  operation  or  during  shut  down. 

B.   INTERNAL  REACTIVITY  FEEDBACK 

In  the  operation  of  a  nuclear  reactor  the  power  generated 
is  directly  proportional  to  the  neutron  population  [1].   In 
order  for  the  neutron  population  to  remain  constant  at  some 
steady  state  power  level,  equation  (2-7)  must  be  equal  to 
zero.   For  this  situation  to  occur,  K,  the  net  internal 
reactivity  which  is  the  sum  of  both  the  reactivity  attributed 
to  the  control  rod  position  and  the  reactivity  due  to  internal 
feedback,  must  be  zero.   Changes  in  the  feedback  reactivity 
depend  upon  temperature  variations  and  fission-product — 
poisoning  which  are  discussed  below 

To  insure  the  stability  of  the  reactor  the  overall  coeffi- 
cient of  reactivity  feedback  must  be  negative  113]  . 
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1.  Temperature  Effects 

When  there  is  a  power  change  in  a  nuclear  reactor 
there  is  considerable  additional  energy  released.   Since  it 
is  not  possible  to  transfer  all  of  the  liberated  energy  to 
the  load,  there  are  temperature  fluctuations  throughout  the 
reactor  and  associated  coolant. 

Temperature  variations  affect  reactivity  in  two 
distant  ways :  first,  the  mean  energy  of  the  thermal  neutrons 
and  hence  their  nuclear  cross-section  vary  with  temperature; 
and,  second,  the  mean  free  path  length  and  leakage  probabilities 
are  functions  of  temperature  [2,  3].   These  effects  result  in 
a  net  negative  reactivity.   A  reactor  power  plant  of  the  type 
described  above,  having  a  negative  temperature  coefficient  of 
reactivity  will  tend  to  be  self-controlling  in  response  to 
load  changes  and  limited  external  reactivity  variations  [14] . 

In  hetrogeneous  reactors  (the  type  under  study  here) 
the  fuel  and  the  coolant  are  at  different  temperatures  and  at 
different  physical  states.   Therefore,  there  will  be  specific 
effects  associated  with  both  the  fuel  and  coolant  [15]. 

Temperature  effects  upon  reactivity  manifest  themselves 
in  temperature  coefficients  of  reactivity,  (a   and  a  ) . 

2 .  Fission-Product  Poisoning 

When  a  reactor  is  operating  at  power  certain  fission 
products,  notably  Xenon  135  and  Samarium  149,  are  accumulated. 
These  products  have  large  nuclear  cross  sections  and  absorb 
the  excess  neutrons,  which  are  vital  in  maintaining  the  chain 
reaction.   The  rate  of  formation  of  these  products  is 
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dependent  upon  the  rate  of  fission  occuring  at  the  specific 
power  level.   Eventually  the  rate  of  formation  and  the  loss 
of  the  absorbing  nuclei  become  equal  and  equilibrium  exists. 
However,  when  the  reactor  is  shut  down,  these  substances, 
which  result  mainly  from  the  decay  of  other  fission  products, 
continue  to  grow  and  may  reach  their  highest  concentration 
several  hours  after  shut  down  occurs .   Since  the  production 
of  these  fission  products  continues  after  shut  down,  it  is 
entirely  possible  that  there  may  be  enough  negative  reactivity 
remaining  to  restrict  the  power  build  up  of  the  reactor  until 
these  products  decay  to  a  level  where  the  inherent  positive 
reactivity  in  the  control  rods  can  overcome  their  adverse 
reactivity  effects. 

The  above  fission  product  phenomenon  does  not  essent- 
ially affect  the  performance  of  a  reactor  operating  at  power. 
For  this  reason  these  effects  are  not  considered  in  the  scope 
of  this  study. 

3 .   Net  Internal  Reactivity 

The  net  reactivity  in  dollar  units,  due  to  both  rod 
position  and  to  temperature  can  be  expressed  as: 

K  =  KR  +  af  TF  +  ac  T^ 

C.   TOTAL  EFFECT  OF  DELAYED  NEUTRONS  AND  NEGATIVE  REACTIVITY 
FEEDBACK 

The  total  effect  of  the  delayed  neutrons  and  negative 

reactivity  can  best  be  illustrated  in  Figure  (3-3) .   These 

plots  are  for  the  specific  plant  under  study.   The  power  plant 
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disturbance  is  a  step  input  of  magnitude  K  =  0.469  dollars 
(i.e. ,  5k  =  0.003)  . 

The  input  causes  the  reactor's  neutron  density  to  rise 
sharply  from  a  steady  state  unit  value  to  a  value  of  1.6.   If 
there  were  no  delayed  neutrons  this  neutron  density  would 
rise  faster  and  obtain  a  larger  value.   The  negative  reactivity 
feedback  acts  rapidly  and  limits  the  excursion  of  the  net 
reactivity  to  a  maximum  value  of  0.38.   After  approximately 
60  seconds,  the  reactivity  feedback  has  adjusted  the  internal 
reactivity  to  zero.   The  steady  state  effect  of  the  disturbance 
is  to  raise  the  neutron  density  to  a  level  of  1.02. 
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IV.   DIGITAL  AND  ANALOG  SIMULATION  OF  POWER  PLANT 
WITHOUT  EXTERNAL  CONTROL 

A.   DIGITAL  SIMULATION 

The  system  equations  were  programmed  and  solved  on  the  IBM 
360-6  7,  using  the  Naval  Postgraduate  School  Library  Program 
(DRKGS)  for  the  integration.   This  integration  method  utilizes 
a  forth-order  "Runge-Kutta"  solution  with  double  precision 
arithmetic. 

Time  response  of  the  system  to  load  changes  of  +40  percent 
are  shown  in  Appendix  B  and  are  tabulated  for  all  experimental 
load  changes  at  the  end  of  this  chapter. 

1.  Load  Changes 

The  system  was  subjected  to  standardized  step  load 
changes  disturbances  of  +  20  percent  and  +  40  percent.   These 
load  changes  were  selected  because  a  20  percent  load  change 
represents  a  change  which  a  maritime  power  plant  should  be 
able  to  respond  to  quickly,  and  a  40  percent  load  change 
represents  a  feasible  change  which  could  be  encountered  during 
an  accidental  situation. 

2 .  Results 

The  results  of  the  simulation  of  this  hypothetical 
system  were  compared  against  the  time  response  of  marine  power 
plant  designs  of  similar  size  [6,  8] .   It  was  found  that  the 
simulation  used  in  this  study  compared  well  with  the  referenced 
systems.   The  steady  state  Cafter  the  load  distrubance)  values 
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of  the  system  variables  coincided  well  with,  their  counterparts, 
The  simulated  system  demonstrated  a  slightly  more  oscillatory 
behavior  than  the  actual  and  proposed  systems  referenced. 
This  oscillatory  behavior  is  attributed  to  the  simplified 
system  model  used  in  this  study.   Even  with  the  small  oscil- 
latory differences,  the  simulation  represents  the  actual 
system  with  a  sufficient  degree  of  accuracy  that  it  is  quite 
suitable  for  an  initial  study  of  marine  reactor  power  plants. 
3 .   Discussion  of  Results 

a)  The  power  plant,  without  benefit  of  an  external 
control  system,  was  shown  to  be  stable  to  all  subjected  load 
changes.   The  system  displayed  a  second  order  well  damped 
response  to  load  changes. 

b)  Temperature  changes,  including  overshoots,  are  not 
severe  and  should  not  cause  thermal  stress  problems. 

c)  System  response,  although  sluggish,  is  fast  enough 
for  use  as  a  maritime  power  plant. 

B.   ANALOG  SOLUTION 

The  analog  simulation  of  the  system  was  carried  out  on 
the  Comcor  CT-5000  computer.   This  computer  has  a  +  100  volt 
dynamic  range.   The  CT-5000  computer  is  interfaced  with  the 
XDS-9300  digital  computer.   This  hybrid  pair  will  be  used 
later  in  the  control  section  of  this  study. 

1.   Prompt- jump  Approximation 

Figure  (B-l)  of  Appendix  B  shows  that  the  time  rate 
of  change  of  N  is  not  large.  For  ease  of  analog  simulation 
(i.e.,  to  reduce  the  number  of  high  gain  amplifiers  required 
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in  the  simulation)  equation  (2-7)  can  be  modified  as  follows 


=  KN  -  N  +  D  (4-1) 


6.4 


making  the  approximation  that  N  =  0  and  solving  for  N 


N  =  JL-  (4-2) 


This  approximation  is  called  the  prompt- jump  approximation 
since  any  change  in  K  is  reflected  instantaneously  in  N.   The 
system  equations  using  this  assumption  were  simulated  on  the 
IBM  360-67  and  it  was  found  that  this  approximation  does  not 
introduce  errors  greater  than  2  percent  in  the  system  response. 
2 .   System  Variables  to  Computer  Variables 

The  system  variables  such  as  temperature  and  relative 
neutron  population  must  be  related  to  analog  computer  variables 
which  are  voltages.   The  system  equations  and  initial  condi- 
tions in  terms  of  computer  variables  are  developed  in  Appendix 
C.   These  equations  are  summarized  below  for  convenience: 
a.   Reactor  Kinetics 

N=frr  (4-3) 

D   =    0.1     (N-D)  (4-4) 


K   =    3.5   KD    -    0.471   T_    -    0.942   T_„  (4-5) 

R  F  AV 
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where 


N(0)   =   10  volts 
5(0)   =   10  volts 


K_(0)  =   35.36  volts 
JK 


b.   Heat  Transfer  From  Fuel  to  Coolant 


TF    =    -    0.5     [Tp    -    TAV]    +    0.5    N  (4-6) 


Where 


TF(0)    =    94.17   volts 


c.   Heat  Balance  for  Coolant 


TAV  =  '5°  TAV.  +  10  TF  +  40  Tc  (4"7) 


Where 


TAV(0)  =  84.17  volts  (4-8) 


d.   Time  Delays 


TBI  "  -°-5  tTBI  "  2  TAV  +  V  (4'9) 


?C  =  °-5  ["fc  +  !*S  +  J*   J  <4-10) 


Where 


TDT(0)  =  86.67  volts 

DX 


T  (0)   =  81.67  volts 
c 
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e.   Heat  Balance  for  Steam  Generator  Secondary-Side 


Tg  =  10.133  TBI  -  0.1324  Tg  -0.0027  \p   Tg]  (4-11) 


Where 

fs(0)  =  79.17  volts 

The  analog  computer  diagram  for  the  solution  of  these  equations 
is  shown  in  Figure  (4-1) . 

3 .   Comparison  of  Analog  and  Digital  Solution 

Below  is  shown  two  tables  comparing  the  results  of  the 
digital  and  analog  computer  solutions.   Because  of  the  inherent 
resolution  quality  characteristics  of  the  digital  computer, 
the  digital  solution  will  be  used  as  the  base  for  determining 
errors.   Settling  time  is  defined  as  the  amount  of  time 
required  for  the  magnitude  of  K(t),  the  net  internal  reactivity, 
to  return  to  a  value  of  +  0.005  after  the  applied  step  load 
disturbances . 

From  the  tabulated  results  it  is  apparent  that  the 
analog  and  digital  solutions  are  in  such  close  agreement  that 
the  analog  solution  is  adaquate  for  the  remainder  of  this  study. 
Henceforth  the  digital  computer  is  used  as  an  auxiliary  compu- 
tational tool  to  solve  specific  problems  but  is  not  used  in 
the  evaluation  of  the  overall  simulation. 

The  value  of  the  system  variables  obtained  from  the 
analog  simulation  were  obtained  through  Analog  to  Digital 
conversion  and  printed  out  from  the  XDS-9300  digital  computer. 
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-12.37 


Fig.  (4-1) . 


ANALOG  COMPUTER  DIAGRAM  FOR  BASIC  POWER 
PLANT. 
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K(MAX/MIN) 
$ 

N(MAX/MIN) 


A  \p    =  +  0 .2 


O 


0.068 


1.21 


En 
M 
O 
M 

Q 


0.066 


1.20     0.83 


A  \i>   =  -0.2 


o 


-0.081 


0.81 


Eh 
H 


-0.084 


0.79 


o 

w 


3.5 


2.53 


N(SS) 


1.21 


1.19 


1.6 


0.82 


0  .80    2.5 


T  (MAX/ 

MIN) 


573.0 


572.7 


0.052 


557.1 


556.8 


0.054 


TF(SS) 


573.0 


572.6 


0.07 


557.4 


557.1 


0.054 


TAV(MAX/ 

MIN) 


499.8 


J500  .5 


0.14 


509.2 


510.1 


0.178 


TAV(SS) 


500.4    501.2 


0.16 


508.2 


509.0 


0.156 


Tg  (MAX/ 


MIN) 


465.3    464.9 


0.086 


486.2 


486.0     0.042 


TS(SS) 

op 


466.0 


.*. 


SETTLING 
TIME 


51.5 


* 


1.199 


465.5 


0.107   I  485.2 
i 


484.9     0.0625 


47.4 


8.6 


1.2 


0.083 


0.801 


•» 


52.0 


0.8 


0.013 


TABLE  (4-1) .   COMPARISON  OF  ANALOG  AND  DIGITAL  SOLUTION 
FOR  A  STEP  LOAD  CHANGE  OF  2  0  PERCENT. 
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A  i>   =   +  0.4 

A  ip  =  -0.4 

* 
o 

3 

3 

Eh 

H 

o 

H 

Q 

O 

o\° 

* 
o 

< 

Eh 
H 
O 
M 

Q 

o 

w 

0\° 

K( MAX/MI N) 
$ 

0.117 

0.119 

1.68 

-0.188 

-  0.19  8 

5.0 

N (MAX /MIN) 

1.40 

1.38 

1.45 

0.586 

0.563 

4.6 

N(SS) 

1.38 

1.37 

0.73 

0.618 

0.595 

4.0 

T  (MAX/ 

MIN) 

580.1 

580.1 

0.0 

548.6 

548.0 

0.123 

Tp(SS) 
oF 

580.1 

579.9 

0.034 

■ 

549.4 

548.8 

0.123 

TAV(MAX/ 

AV    MIN) 
oF 

495.7 

496.3 

0.12 
„ ,,.„,, 

515.0 

516.2 

0.232 

TAV(SS) 

op 

497.0 

497.5 

0.10 

512.2 

513.0 

0.156 

T  (MAX/ 

MIN) 
oF 

456.0 

455.3 

0.154    498.0 

495.0 

0.60 

Tg(SS) 
oF 

457.0 

456.4 

0.133 

495.2 

495.2 

0.0 

SETTLING 
TIME 

54.2 

50.0 

8.4 

108.0 

98.0 

10.0 

* 

1.397 

1.4 

0.215 

0.604 

0.6 

0.667 

TABLE  (4-2) .   COMPARISON  OF  ANALOG  AND  DIGITAL  SOLUTION  FOR 
A  STEP  LOAD  CHANGE  OF  40  PERCENT. 
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V.   CONTROL 

The  purpose  of  an  external  control  scheme  for  a  nuclear 
power  plant  must  be  to  set  up  the  desired  steady  state 
operating  conditions,  and  to  restrict  transients  in  vital 
plant  variables  to  be  consistant  with  design  limitations. 
The  above  purposes  should  be  obtained  without  any  significant 
reduction  in  the  inherent  stability  characteristics  of  the 
plant  without  external  control. 

Examination  of  the  basic  plant  transient  responses  in 
Appendix  B  indicates  that  the  system  under  study  is  very 
stable  when  subjected  to  load  changes.   It  is  however, 
desirable,  as  pointed  out  below,  to  limit  the  excursion  of 
certain  plant  variables  from  their  steady  state  normal  value. 
In  a  plant  with  a  large  degree  of  inherent  stability  direct 
control  of  plant  variables  may  be  attempted.   If  a  plant  has 
poor  inherent  stability,  neutron-level  control  should  also  be 
used  [3] . 

A.   STEADY  STATE  PROGRAMMING 

The  pattern  that  the  plant  temperatures,  pressures,  and 
flow  rates  assume  as  a  function  of  power-demand  is  oftened 
referred  to  as  a  program  [4] .   A  specific  control  strategy 
can  be  devised  to  obtain  the  desired  program.   The  two  most 
significant  programs  for  discussion  are  the  constant  T  „  and 
constant  T_  programs. 
/ 
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1.      Constant  T 


AV 


TEMPERATURE 
OR 
PRESSURE 


POWER  DEMAND 

Fig.     (5-1) .    SKETCH    OF    TEMPERATURE    AND    PRESSURE 
PROFILES    FOR   A   CONSTANT    T         PROGRAM. 

nV 

Figure  (5-1)  is  a  graphical  representation  of  the 
constant  T    program.   At  zero  power  demand,  all  significant 
plant  temperatures  (except  T^)  are  at  the  same  value.   As  the 

r 

power  demand  increases,  the  steam  temperature  and  hence 
pressure  decrease,  while  T   and  T   vary  in  such  a  manner  as 
to  keep  T    constant. 


2.   Constant  T 


TEMPERATURE 
OR 
PRESSURE 


POWER  DEMAND 


Fig.     (5-2) .    SKETCH   OF    TEMPERATURE    AND    PRESSURE 
PROFILES    FOR   A   CONSTANT    T      PROGRAM. 
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Again  at  no  load  all  pertinent  temperatures  are  equal. 
As  the  load  increases  the  temperature  and  the  pressures  follow 
the  pattern  as  shown  in  Figure  C5-2) . 

3.   Comparison  of  Constant  T  v  and  Constant  T   Programs 

The  constant  T  „  program  is  the  one  preferred  by  the 
reactor.   The  principle  advantage  of  the  constant  T   program 
for  this  type  reactor  is  that  the  primary  coolant  pressure 
also  remains  constant.   Therefore  the  need  for  an  external 
pressurizer  in  the  primary  circuit  is  reduced. 

One  additional  advantage  occurs  if  the  plant  under 
study  had  been  fuel  with  an  uranium  alloy  vice  an  oxide  of 
uranium.   For  a  metal  alloy  fuel  the  feedback  of  internal 
reactivity  associated  with  changes  in  fuel  temperature  is 
insignificant  when  compared  to  the  feedback  associated  with 
changes  in  average  coolant  and  therefore  can  be  ignored  [3] . 
For  this  case  equation  (2-9)  becomes: 

K  ■  KR  +  °c  TAV  f5"1* 

Since  in  stead  stake  K  must  equal  zero,  it  is  apparent 
from  equation  (5-1)  that  there  is  no  need  for  control  rod 
motion  to  maintain  T  constant.  However  the  plant  under  study 
is  fueled  with  UO-  and  changes  in  fuel  temperature  have  a 
strong  affect  on  net  reactivity.  Examination  of  equation 
(2-9)  now  indicates  that  there  is  a  requirement  to  change 
control  rod  position  to  keep  the  average  temperature  constant. 

The  advantages  of  the  constant  Tq  Program  are  apparent, 
from  Figure  (5-2) ,  in  the  design  of  the  secondary  loop  of  the 
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plant.   A  constant  Tg  program  also  provides  a  constant  steam 
pressure  which  allows  for  an  optimum  design  of  the  steam 
plant  with  conventional  controlling  devices. 

The  two  primary  disadvantages  of  the  constant  Tq 
program  are  first,  allowing  T   to  vary  over  a  large  range 
(required  to  keep  T„  constant)  may  cause  serious  problems  in 
the  primary  coolant  pressure  system  and  second,  when  the  plant 
is  changing  from  one  steady  state  power  level  to  another,  the 
plant  must  fight  the  tendency  of  the  coolant  negative  temper- 
ature coefficient  to  hold  T,TT  constant. 

AV 

With  the  exception  of  Nautilus,  which  uses  a  constant 
Tp  program,  all  Naval  plants  use  a  constant  T  v  program  [16] . 
The  N.  S.  Savannah  utilizes  a  control  plan  for  a  constant  T __. 
program  with  a  plus  or  minus  3°F  dead  zone  about  the  normal 
load  T    value  [17] . 

B.   DESIGN  OF  THE  EXTERNAL  CONTROL  SYSTEM 

Consistant  with  the  advantages  of  the  constant  T    program, 
the  simulated  power  plant  will  be  controlled  under  a  constant 
TAV  program.   As  noted  in  section  (5A-3)  control  rod  movement 
is  required  in  order  to  obtain  the  desired  program.   The  next 
decision  to  be  made  is  at  what  rate  may  reactivity  be  externally 
introduced  into  the  reactor.   The  answer  to  this  and  other 
design  questions  are  discussed  below. 

1.   Reactivity  Rate 

The  one  basic  underlying  consideration  which  dictates 
the  design  of  a  control  system  for  a  nuclear  power  plant  is 
safety.   As  shown  in  Figure  (3-3)  the  power  plant  operating 
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at  normal  power  can  withstand  a  step  reactivity  input  of 
magnitude  5k  =  .469$.   This  however  is  not  the  case  for  a 
reactor  operating  at  reduced  power  levels  where  the  negative 
reactivity  feedback  effect  is  reduced.   Therefore  in  order  to 
provide  a  maximum  safety  margin  to  avoid  startup  or  low  power 
accidents  input  reactivity  rates  will  be  restricted  to  a  value 
which  can  be  handled  by  a  reactor  operating  at  lower  power. 
The  solution  to  equations  (2-7)  and  (2-8)  for  ramp  input  of 
reactivity,   K(t),  does  not  give  a  valid  solution  to  the 
problem  because  there  is  no  feedback  considered.   Externally 
introduced  reactivity  rates  are  determined  largely  from 
empirical  data  on  the  individual  plants  involved.   A  reasonable 
reactivity  rate  for  a  plant  of  this  size  is  IOC/sec  (0.1  $/sec) 
[71. 

It  is  feasible  that  this  rate  of  reactivity  may  be 
controlled  by  a  manual  means.   However,  it  was  desired  that 
the  system  be  fully  automated  therefore  an  automatic  controller 
was  used. 

2 .   Automatic  Controller 

The  automatic  controller  consists  of  a  hybrid  combina- 
tion of  an  analog  control  rod  positioning  system  and  a  digital 
system  for  comparison  of  T    plus  the  time  derivative  of  T  v 
with  the  reference  value  of  505°F.   It  was  decided  to  use  the 
+  3°F  dead  band  error  detector  of  N.  S.  Savannah  as  a  valid 
method  to  obtain  a  discontinuous  controller.   The  principle 
advantage  of  a  discontinuous  controller  is  that  it  utilizes 
the  reactor  inherent  stability  characteristics  to  reduce  wear 
on  the  control  rod  mechanisms.   The  dead  zone,  although  it 
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must  be  consistent  with,  the  design  temperature  (.i.e.,  pressure) 
fluctuations  in  the  primary  loop,  discriminates  against  noise 
on  the  T    signal.   A  block,  diagram  description  of  the  entire 
simulation,  plant  plus  controller,  is  represented  in  Figure 
(5-3) .   The  design  of  the  control  rod  positioning  system  and 
a  discussion  of  the  hybrid  simulation  along  with  transient 
response  curves  for  a  step  load  input  of  +40  percent  of  the 
system  with  external  control  are  shown  in  Appendix  D. 
3.   Evaluation  of  External  Controller 

The  external  controller  was  implemented  and  incorpor- 
ated into  the  simulation.   The  system  was  subjected  to  step 
load  disturbances  of  plus  and  minus  20  and  40  percent.   The 
improvement  in  system  performance  can  be  seen  from  tables  (5-1) 
and  (5-2) .   The  reduction  of  the  excursions  of  T   ,  both 
transient  and  steady  state,  and  the  reduction  of  settling  time 
were  considered  significant. 

The  plant  was  subjected  to  step  load  changes  of  up  to 
80  percent  (an  unrealistic  change) .   It  was  found  that  the 
system  was  stable  and  that  the  system  response  improved  up 
to  step  load  inputs  of  60  percent.   Above  60  percent,  the 
system  response  deteriorated  and  eventually  went  into  a  limit 
cycle  when  the  step  load  disturbance  became  greater  than  75 
percent.   The  deteriorated  performance  and  eventual  instability 
of  the  system  were  attributed  to  the  fact  that  the  derivative 
of  T  v  became  very  large  and  began  to  work  in  the  opposite 
direction  of  the  T  v  excursions  rather  than  in  normal  anti- 
cipatory mode  of  a  position  plus  derivative  feedback  controller 
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A    4>    =    +0.2 

i    f   =    -0.2 

EXTERNAL 
CONTROLLER 

NO 

EXTERNAL 
CONTROLLER 

EXTERNAL 
CONTROLLEF 

NO 

EXTERNAL 

CONTROLLER 

MAV 

-      3       °F 

-       5.2°F 

+    4.2°F 

+       4.2°F 

TAV(SS) 

505. 6°F 

500.  4°F 

505       °F 

508. 2°F 

K(MAX/MIN) 

0.12$ 

0.068 

-    0.16$ 

-      0.081 

AKR 

+      0.1$ 

0.0 

SETTLING 
TIME 

36    sec 

51.5    sec 

40.0    sec 

49.5    sec 

TABLE  (5-1)  .   COMPARISON  OF  SIMULATED  POV7ER  PLANT 
WITH  AND  WITHOUT  EXTERNAL  CONTROLLER 
FOR  A  STEP  LOAD  CHANGE  OF  20  PERCENT. 
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a  i>  = 

+  0.4 

a  i>  =  - 

•0.4 

EXTERNAL 
CONTROLLER 

iv  i — i 

oSO 

2  W  h 

So 

EXTERNAL 
CONTROL  LEF 

NO 

EXTERNAL 
CONTROLLEF 

AT 
AV 

-      5.1 

-    9.3°F 

+       7.2°F 

+10.0°F 

TAV(SS) 

508. 0°F 

497. 0°F 

507. 4°F 

512. 2°F 

K( MAX/MI N) 

0.17$ 

0.117$ 

-      0.32 

-    0.188$ 

AKR 

+      0.6$ 

+      0.6$ 

SETTLING 
TIME 

36    sec 

54.2    sec 

77    sec 

108    sec 

TABLE  (5-2) 


COMPARISON  OF  SIMULATED  POWER  PLANT 
WITH  AND  WITHOUT  EXTERNAL  CONTROLLER 
FOR  A  STEP  LOAD  CHANGE  OF  40  PERCENT, 
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The  controller  was  considered  adequate  for  all 
reasonable  power  changes. 
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VI .   CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  and  recommendations  are  offered 
as  a  result  of  the  investigation  contained  in  this  report: 

A.   CONCLUSIONS 

1.  A  properly  designed  pressurized  water  nuclear  reactor 
(i.e.,   one  with  a  negative  temperature  coefficient  of 
reactivity)  possesses   a  great  deal  of  inherent  stability  when 
the  associated  power  plant  is  operating  in  the  normal  power 
range.   It  was  observed  that  the  negative  temperature  coeffi- 
cient of  reactivity  served  to  stabilize  the  reactor  for  dis- 
turbances in  both  load  and  input  reactivity. 

2.  The  naturally  occuring  phenomenon  of  delayed  neutrons 
greatly  simplifies  the  control  problems  associated  with 
nuclear  reactors.   The  effect  of  this  phenomenon  is  more 
evident  at  the  lower  power  levels  before  the  temperature 
reactivity  feedback  effect  has  become  affective. 

3.  The  control  of  a  stable  nuclear  power  plant  operating 
in  the  power  range  under  a  constant  T    program  with  constant 
coolant  flow  rate  is  a  relatively  simple  task.   If  sufficient 
stability  exists  in  the  reactor  then  direct  control  of  a  plant 
parameter  (T   )  can  be  attempted. 

4.  If  step  load  changes  greater  than  60  percent,  for  which 
the  above  control  system  is  not  suitable,  are  considered 
feasible  a  comparator  could  be  installed  in  the  differentiating 
circuit  to  disregard  the  derivative  if  it  becomes  too  large. 
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B.   RECOMMENDATIONS 

Further  workers  in  this  area  should  consider  the  following 
facets  of  reactor  systems. 

1.  The  above  study  was  concerned  with  the  reactor  power 
plant  dynamics  when  excited  by  both  load  disturbances  and 
changes  in  reactivity  investigated  separately.   Another  area 
of  interest  is  the  response  of  the  reactor  under  the  condition 
of  load  changes  and  externally  induced  reactivity  changes 
(both  aiding  and  opposing)  occuring  simultaneously.   This 
situation  might  represent  an  accidental  condition. 

2.  The  T  v  program  used  to  control  the  power  plant  does 

T„   +  T 
H     c 

not  consider  the  maximum  value  to  TTT.   Since  T,TT  =  ~ 

H  AV      2 

the  value  of  T„  might  become  too  high  for  design  considera- 
tions.  To  institute  a  control  program  with  a  constant  T  v 
and  also  limiting  T„  to  some  maximum  value  would  make  the 
requirement  of  a  variable  coolant  flow  rate  necessary.   This 
added  requirement  would  necessitate  a  new  mathematical  model. 
Under  the  assumption  of  a  constant  flow  rate,  the  time  coolant 
spends  in  the  reactor  (xn)  is  a  constant.   However,  in  the 
model  for  a  plant  where  a  maximum  value  of  T„  must  be 


H 


considered  t  ~   would  be  a  function  of  T„ . 

U  n 


3.   More  sophisticated  control  problems,  than  reactor 
control  at  operating  conditions,  exist  in  the  nuclear  power 
plant  control  field.   These  are  the  problems  associated  with 
the  reactor  in  the  startup  range  and  problems  involved  in  the 
mechanical  control  of  the  various  plant  components  of  the 
secondary  loop.   The  problems  in  the  control  of  the  secondary 
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loop  result  from  the  fact  that  when  the  reactor  is  controlled 
under  a  constant  T_„  program,  the  steam  temperature  and 
pressure  vary  significantly  over  the  operating  range  of  the 
plant.   It  is  recommended  that  further  investigators  consider 
the  control  problems  of  the  secondary  loop. 
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APPENDIX  A 
SPECIFIC  SYSTEM  EQUATIONS 

In  the  development  of  the  generalized  system  equation, 
equations  (2-28)  and  (2-34)  are  expressed  in  terms  of  the 
variable  T_. .   T_.  is  not  a  state  variable  of  the  system.   The 
equation  for  the  primary  side  of  the  heat  exchanger  (2-31)  can 
be  solved  for  TR  and  this  value  substiuted  into  the  equations 
(2-28)  and  (2-34)  to  place  them  in  a  more  suitable  form. 
This  yields : 

•  /  2T 

m  _      _       1 |         rp  _      S_  E~l 

c  t,    I       C        e+1        E+l      BIj 


4 


Ts  =  -M    "     eTT      tbi  +    [£i  +  kt  *  Sf 


)■■) 


The  system  parameters  for  the  hypothetical  plant  are  stated 
below.   It  should  be  emphasized  that  these  parameters  do  not 
represent  any  actual  real  system  but  that  the  parameters  used 
here  are  obtained  from  various  sources  in  the  literature  [5, 
6,  7 ,  8,  9]  and  represent  typical  values  which  might  be 
associated  with  a  plant  the  size  of  the  power  plant  of  N.  S. 
Savannah. 

1)  Fraction  of  delayed  neutrons    3  =  0.0064 

2)  Average  delay  constant  of  the  6  group  delayed 
neutrons    X  =  0 . 1 
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-3 

3)  Mean  neutron  life  time    %   =  10    sec 

4)  Steady  state  temperature 

A)  Average  fuel  temperature      TF(0)  =  56  5 °F 

B)  Average  coolant  temperature   T  v(0)  =  505°F 

C)  Boiler  inlet  temperature     TDX(0)  =  520°F 

D)  Steam  temperature  Tg  (0)  =  475°F 

E)  Cold  leg  temperature  Tc^  =  490°F 

5)  Temperature  coefficient  of  reactivity  for  the  fuel 

(uo2) 

ciF  =  -1.57  x  10"*2  $/°F 

6)  Temperature  coefficient  of  reactivity  for  the  coolant 
(pressurized  water) 

a   =  -3.14  x  10~2  $/°F 
c 

7)  Time  constants  for  the  various  heat  transfer  equations 

A)  T,  =   2.0  sec. 

B)  t2  =   0.1  sec. 

C)  t_  =   2 .0  sec* 

D)  t.  =   2.0  sec* 

E)  tc  =   5.0  sec* 

8)  In  the  region  of  interest  the  steam  pressure-steam 
temperature  relation  is  approximately  constant 

2PS 

=  5  psi/°F 


2T 
Z1S 


9)   The  values  of  K  (0) ,  A,,  KT/  x  ,  t 2  were  determined 

by  setting  the  equation  in  which  these  terms  appear  to 
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zero  (steady  state)  and  solving  for  the  unknown 
variables.   The  results  are: 
A)   KR(0)  =  24.7275  $ 


B) 

Al 

= 

30.0 

C) 

E 

= 

2.0 

D) 

KT 

= 

1.33   x    10"2 

E) 

To 

= 

0.05    sec. 

F) 

To 

= 

0.02    sec 

*T3'  T4'  an<^  T5  were  obtained  by  engineering  approxi- 
mation and  are  substantially  less  than  those  mentioned 
for  a  stationary  power  plant  [9]. 
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APPENDIX  B 
TIME  RESPONSE  OF  POWER  PLANT 

As  a  representative  example  of  system  time  response 
Figures  (B-l) ,  (B-la) ,  (B-2)  and  (B-2a)  of  this  appendix  show 
the  time  response  of  selected  system  variables  to  a  step  load 
change  of  +40  percent  normal  load.   Figures  (B-l)  and  (B-la) 
are  solutions  using  the  digital  simulation.   Figures  (B-2) 
and  (B-2a)  are  solutions  utilizing  the  analog  simulation. 
Pertinent  values  of  the  variables  are  shown  adjacent  to  their 
time  response.   The  system  response  for  other  experimental 
load  changes  are  similar  in  nature  to  those  shown. 
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Fig.  (B-l) 


TIME  RESPONSE  OF  SIMULATED  POWER 
PLANT  TO  A  40  PERCENT  STEP  LOAD  CHANGE 
(DIGITAL) . 
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T  (MIN)  =  455. 3 °F 
TS(SS)   =  456.  4°F 


0         40         80         120 

STEAM  TEMPERATURE  VERSUS  TIME 
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0        40         80         120        160  sees 

AVERAGE  COOLANT  TEMPERATURE  VERSUS  TIME 


Fig.  (B-la) .   TIME  RESPONSE  OF  SIMULATED  POWER  PLANT 

TO  A  40  PERCENT  STEP  LOAD  CHANGE  (DIGITAL) 
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TIME  RESPONSE  OF  SIMULATED  POWER  PLANT 
TO  A  40  PERCENT  STEP  LOAD  CHANGE  (ANALOG) 
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Fig.  (B-2a) .   TIME  RESPONSE  OF  SIMULATED  POWER  PLANT 

TO  A  40  PERCENT  STEP  LOAD  CHANGE  (ANALOG) 
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APPENDIX  C 

PREPARATION  OF  SYSTEM  EQUATIONS 
FOR  USE  ON  ANALOG  COMPUTER 


Equations  (4-2)  and  (2-36)  through  (2-42)  describe  the 
power  plant  in  terms  of  system  variables.   These  system  vari- 
ables have  units  such  as  degrees  Fahrenheit  and  dollar  units 
of  reactivity  which  must  be  related  to  analog  computer  vari- 
ables which  are  voltages.   In  addition  to  the  conversion  of 
units,  the  analog  voltages  have  a  maximum  value  of  100  volts 
so  scaling  of  the  variables  may  also  be  required. 

To  accomplish  this  conversion  a  scaling  factor  coefficient 
x  is  introduced.   The  relation  between  computer  variables  and 
system  variables  then  becomes: 


N   = 


Vs 


D   =   XDD 


T    =   v    T 
XF      *TF   F 


and  so  forth. 

The  individual  x's  are  equal  to  the  maximum  expected  value 
of  the  system  variable  divided  by  the  maximum  desired  voltage 
for  the  computer  variables.   For  example: 

TF   =   XTF  TF 
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assuming  T_  (MAX)  =  600°F  and  T_  (MAX)  =  100  volts 

TF 
-   x   =   =£-  =  6°F/volt 

F 
A  summary  of  the  scaling  coefficients  is  as  follows: 

xN  =  .1  volt 

x_  =  .  1  volt 
xK  =  .2$/volt 

XKR  =  *7  $/volt 
X  ip  =  .  ]  volt 

XTF  =  XTS  =  ^AV  =  *TC  =  XTBI  =  6°F/volt 

The  system  equations  are  now  altered  by  replacing  the  system 

variables  with  their  corresponding  computer  variables  and 

conversion  factors.   The  initial  condition  of  the  differential 

equations  are  determined  as  follows : 

_       T  (0) 
T^(0)  = 


F        *TF 


In  a  final  form  the  equations  now  become 
1.   Reactor  Kinetics 


i-  52. 


5-K 

D  =  0.1  (N-D) 


K  =  3.5  K   -  0.471  TF  -  0.942  TAy 
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Where 

N(0)  =  10  volts 
5(0)  =  10  volts 
KR(0)  =  35.36  volts 
2.   Heat  Transfer  from  Fuel  to  Coolant 


Tp  =  -  0.5  (Tp  -  TAV)  5.0  N 


Where  TF(0)  =  9  4.17  volts 
3.   Heat  Balance  for  Coolant 


T    =  -50  T    +  T   +  40  T 
AV     DU      AV     F    H      XC 


Where  TAV(0)  =  8  4.17 


4.   Time  Delays 


TBI  =  °-5  ("TBI  +  2TAV  "  V 


Tc  =  0.5  (-Tc  +  2/3  Tg  +  1/3  TBJ) 


Where  TDT(0)  =  86.67  volts 


Tc(0)  =  81.67  volts 


5.   Heat  Balance  for  Steam  Generator  Secondary-Side 
Ta    =  (0.133  tTT  -  .1324  tI  -  0.0027  T~J 


Where  Tg(0)  =  79.17  volts 


$(0)  =   5  volts 
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APPENDIX  D 
CONTROLLER  SIMULATION 

Analytical  design  of  the  external  controller  could  have 
been  approached  in  one  of  two  ways.   First  a  small  signal 
linearization  of  the  plant  dynamics  could  have  been  attempted. 
This  method  was  not  valid  since  two  of  the  plant  variables, 
N  and  D,  change  in  almost  direct  proportion  to  load  changes. 
For  example  a  40  percent  load  change  causes  a  40=  percent 
change  in  N  and  D.   This  change  in  magnitude  of  the  variables 
was  outside  the  range  of  small  signal  application.   Second  a 
describing  function  could  have  been  obtained  utilizing  sinu- 
sodial  inputs  [18].   The  describing  function  technique  is  a 
laborious  method  of  obtaining   insight  into  the  absolute 
stability  limits  of  the  system  but  gives  little  insight  into 
the  controller  design  problem.   In  this  initial  study  the 
controller  design  and  evaluation  was  carried  out  through 
simulation.   If  a  more  detailed  study  had  been  undertaken,  a 
describing  function  would  have  been  necessary  to  verify  the 
simulation  results. 

The  simulation  of  the  external  controller  was  accomplished 
by  use  of  a  trial  and  error  technique  utilizing  the  hybrid 
computing  feature  of  the  XDS-9  300  and  Comcor  CI-500  0  computers, 

The  analog  voltage  representing  T    was  sampled  at  a 
frequency  of  10  samples  per  second  (this  sampling  frequency  is 
well  in  excess  of  that  required  by  the  sampling  theorem) .   The 
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values  of  T   were  stored  in  an  array.   When  the  above  array 
had  accumulated  10  members/  an  approximation  to  the  time 
derivation  of  T   (10)  -  T   (1)  (divided  by  unit  time) .   The 
inclusion  of  the  derivative  feedback  will  give  the  system  an 
anticipatory  response.   After  the  subtraction  the  array  was 
updated  so  that  a  new  derivative  was  available  every  0.1 
second.   Due  to  the  fact  that  TAV  was  changing  at  a  rate  of 
less  than  0.5  degrees  per  second  at  distrubances  below  60 
percent,  the  approximation  of  the  derivative  of  T  v  was  con- 
sidered  accurate.   The  sum  of  TAV  and  T  v  were  subtracted  from 
T  v  references  (505°F) .   If  the  error  were  more  than  3°F  an 
error  signal  of  +  1  volt  was  sent  to  the  analog  computer. 


+  V 


K 


m 


S(St  +1) 
m 


+  0 


WORM  GEAR 
ASSEMBLY 


±  &Kn 


Fig.  (D-l) .   BLOCK  DIAGRAM  OF  CONTROL  ROD 
POSITIONING  SYSTEM. 


The  +  1  volt  error  signal  was  placed  into  a  standard 


K 


electrical  motor  with  transfer  function 


m 


S(St   +  1)  * 
n 


As 


indicated  in  Figure  (D-l)  in  the  actual  plant  this  motor  would 
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drive  at  constant  speed  through  a  worm  gear  assembly  and 
remove  or  insert  the  control  rods  at  a  constant  speed.   The 
gear  assembly  represents  only  a  linear  gain  and  can  be  com- 
bined with  the  motor  gain  to  give  a  new  transfer  function 

q  #e    ,  -■  \  •   The  control  rod  effectiveness  (reactivity  per 
m 

inch)  is  a  non-linear  function  of  position  and  is  different 

for  individual  reactors.   It  is  assumed  here  that  the  rods 

are  operating  in  a  linear  mode  where  <5K  =  .1  $/sec.   Under 

the  assumption  that  x   =  .1  sec.  the  final  transfer  function 

m 

becomes  g  ,g  +  T q\     •   Simulation  of  the  drive  motor  yields  the 

transient  responses  of  Figures  (D-2)  and  (D-3) .  From  the 
simulation  of  the  control  rod  positioning  device  the  rise 
time  is  0.25  seconds  and  the  steady  state  6K   is  10.15^/sec. 

The  rise  time  could  be  shortened  and  the  steady  state  K 
remain  constant  if  the  forward  gain  were  increased  and  a 
tachometer  feedback  loop  used.   However  this  would  allow  the 
possibility  of  an  excessive  6K   if  the  tachometer  loop  were  to 
fail.   Since  the  rise  time  of  the  uncompensated  system  was 
quite  satisfactory  for  this  application  the  tachometer  feedback 
loop  was  not  employed. 

Transient  response  of  the  system  with  external  controller 
is  shown  in  Figure  (D-4)  and  (D-4a) . 
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